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Abstract. This paper presents results concerning bifurcations of 2D 
piecewise— smooth dynamical systems governed by vector fields. Generic 
three— parameter families of a class of Non— Smooth Vector Fields are 
studied and the bifurcation diagrams are exhibited. Our main results 
describe the unfolding of the so called Resonant Fold — Saddle singu- 
larity. 



1. Introduction 

Bifurcation theory describes how continuous variations of parameter va- 
lues in a dynamical system can, through topological changes, cause the 
phase portrait to change suddenly. In this paper we focus on certain unsta- 
ble non— smooth vector fields within a generic context. The framework in 
which we shall pursue these unstable systems is sometimes called generic 
bifurcation theory. In [1] the concept of A;*'^- order structural stability is 
presented; in a local approach such setting gives rise to the notion of a 
codimension k singularity. Observe that, so far, bifurcation and normal form 
theories for non— smooth vector fields have not been extensively studied in 
a systematic way. 

Non— smooth dynamical systems (abbreviated by NSDS) has become cer- 
tainly one of the common frontiers between Mathematics and Physics or 
Engineering. Problems involving impact or friction are piecewise— smooth, 
as are many control systems with thresholds. In this article the bifurcation 
diagrams of some typical singularities of NSDS in the plane are discussed. 
We study in this setting a set of typical bifurcations which are not found 
in smooth systems. We focus our attention on Filippov systems (see [7]), 
which are systems modelled by ordinary differential equations with discon- 
tinuous righthand sides. It is well known that many of these models (see 
for instance [2] and [3]) occur in generic fc— parameter families and therefore 
they typically undergo generic codimension k bifurcations. Many authors 
have contributed to the study of Filippov systems (see for instance [7| and 
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|10j). One of the starting points for a systematic approach in the geomet- 
ric and quahtative analysis of NSDS is the work [13], of M. A Teixeira, on 
smooth systems in 2— dimensional manifolds with boundary. The generic 
singularities that appear in NSDS, as far as we know, were first studied in 
|14j . Bifurcations and related problems involving or not sliding regions are 
studied in papers like [5] and [8]. The classification of codimension 1 local 
and some global bifurcations for planar systems was given in [TT]. In [9] 
is shown how to construct the homeomorphisms which lead to equivalences 
between two non— smooth systems when the discontinuity set is a planar 
smooth curve. In that work codimension two singularities were discussed 
and amazing phenomena in their bifurcation diagrams appeared in the form 
of infinitely many branches of codimension 1 global bifurcations. These bi- 
furcations, that also appear in the present work, called ST— bifurcations, 
are characterized by the connection between a saddle critical point and a 
tangency singularity. See [15] for a survey on NSDS and references there in. 
Those papers give the necessary basis for the development of our approach. 

The specific topic addressed in this paper is the characterization of spe- 
cific families of the Resonant Fold— Saddle bifurcation diagram. 

Let C be a compact set and T, K given by S = f~^{0), where / : 
X — 7- M is a smooth function having G M as a regular value (i.e. 'Vf{p) ^ 0, 
for any p € /~^(0)) such that dKCiT, = or dK rti S, where dK is a smooth 
manifold. Clearly the switching manifold S is the separating boundary of 
the regions S+ = {g G K\f{q) > 0} and S_ = {g G K\f{q) < 0}. We 
can assume that S is represented, locally around a point q = {x,y), by the 
function f{x,y) = y. 

Designate by the space of C— vector fields on K endowed with the 
C^— topology with r > 1 or r = oo, large enough for our purposes. Call 
= Q^[K, f) the space of vector fields Z : K such that 



where X = {fi,gi), Y = (/2,52) are in x^- We write Z = {X,Y), which 
we will accept to be multivalued in points of S. The trajectories of Z are 
solutions of = Z{q), which has, in general, discontinuous right— hand side. 
The basic results of differential equations, in this context, were stated by 
Filippov in [7] . Related theories can be found in [lOl [I2l [14] . 

Definition 1. Two non— smooth vector fields Z, Z £ i}^{K,f) defined in 
open sets U, U £ K and with switching manifolds Ti C U and Ti C U 
respectively are T,— equivalent if there exists an orientation preserving ho- 
meomorphism h : U ^ U which sends T, to T, and sends orbits of Z to orbits 




for (x,?/)gS+, 
for (x,y)GS_, 



ofZ. 
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We say that two unfoldings Ba : x M'^ ^ ^^j^^ : x ^ ]R2^ 
where A G and fi G M',0, are topologicahy equivalent if for each 

A G M*^,0 there exists A{X) G R\0 such that the vector fields 9a and 
^A{A) ^re S— equivalent. And we say that an unfolding 0a is generic if in a 
neighborhood of ©a any other unfolding is topologicahy equivalent to ©a- 

In what follows we will use the notation 

X.fip) = {\/f{p),X{p)) and Y.fip) = {V f{p),Y{p)) . 

1.1. Setting the problem. Let Xq be a smooth vector field defined in S+. 
We say that a point po G 5] is a T,—fold point of Xq if Xo./(po) = but 
Xq./(po) 7^ 0. Moreover, po £ 5] is a visible (respectively invisible) S— fold 
point of Xo if Xo.fipo) = and Xlf{po) > (resp. Xlf{po) < 0). We 
denote the set of all vector fields defined in S+ presenting a S— fold point by 
r^^. We endow with the C^— topology. In this universe a S— fold point 
has codimension zero. It is possible to consider /(x, y) = y and the following 
generic normal forms Xo(x,y) = (ai,/3ix) with ai = ±1 and (3i = ±1 (see 
|16| . Theorem 2). 

Let Yo be a smooth vector field defined in S_. Assume that Yq has a 
hyperbolic saddle point Sy^ on S and that the eigenspaces of DYq{Syo) 
transverse to S at Syq- We denote the set of all vector fields defined in S_ 
presenting a hyperbolic saddle with the eigenspaces transverse to S by . 
We endow Tj^_ with the C""— topology. In this universe a saddle point Syq 
has codimension one. We say that two vector fields Y,Y ^ Fj, defined in 

open sets U and U, respectively, are C'^— orbitally equivalent if there exists 
an orientation preserving homeomorphism h : U ^ U that sends orbits of 
Y to orbits of Y. Prom jl3j we know that any saddle 5*0 is orbitally 
equivalent to its linear part by a S— preserving homeomorphism. And the 
linear saddle with eigenspaces transverse to the x— axis has the generic nor- 
mal forms Yo{x,y) = {a2y,a2x) with 02 = il- So the generic unfolding of 
the singularity is given by Yg = (02 (2/ + P),02x) where /3 G M. 
Let U he a small neighborhood of Yq in F^ • Then: 

(a) There exists a smooth function L : [/ — t- M, such that DLy^ is surjective. 

(b) The correspondence Y — >• Sy is smooth, where Sy is a saddle point of 
Y. 

(c) If L(y) > then 5y G S_. 

(d) If L{Y) = then Sy G S. 

(e) If L(y) < then 5y G S+. 

In this paper we are concerned with the bifurcation diagram of systems 
Zq = (Xq, Yq) in QJ' such that pQ = Sy^ G S. This singularity will be called 
Fold — Saddle singularity (see Figures [T] and [2] — the dotted lines in these 
and later figures represent the points where X.f = and Y.f = 0). 
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Figure 1. (Invisible) 
Fold- Saddle Singularity. 



Figure 2. (Visible) 
Fold-Saddle Singularity. 



Let p = (0,0) be a fold— saddle singularity of Z = {X,Y). We denote 
the set of all non— smooth vector fields Z = {X, Y) such that X G F^^ 

by F^~'^. We endow F^"*^ with the product topology. Let 



and G F| 



Zo = {Xo,Yo) e F^"^. Observe that is the unique singularity of Xq 
around a neighborhood Wo of the origin in M^. There exists a neighborhood 
Uq of Zq in such that for any Z = (X, Y) £ Uq we may find a S— fold 
point pz = ikz,0) G Wo such that it is the unique singularity of X in Wq. 
Moreover the correspondence Z ^ pz is C"'. 

In the same way, for any Z = {X, Y) £ Uq we find a C^— correspondence 
B : Uq ^ where B{Z) = sz = iaz,bz) is the (unique) equilibrium 
(saddle) of Y in Uq. We are assuming that the eigenspaces of DYs^{qz) are 
transverse to S at sz- We have to distinguish the cases: (i) bz < 0, [ii] 
bz = and {iii) bz > 0. Observe that when bz < (resp. bz > 0) there 
is associated to Z an invisible (resp. visible) S— fold point of Y given by 
Qy = (cz,0) e Wq. Moreover lirrn,^^QCz = az- 
Define F{Z) = {kz - az, bz). We get: 



(i) The derivative DF : Uq 



p2 ; 



IS surjective; 



(ii) F ^(0) = 0,2 is a codimension two submanifold of rj*". 

Then this fold— saddle singularity occurs generically in two— parameter 
families of vector fields in . 

We consider the following model: 



(1) 



X' 



Y 



±1 

ai{T)x 

hy 

kix 



if y > 0, 
if y < 0, 



where ai{inv) = —1, ai{vis) = 1 and ki = ±1. 
Lemma 1. If Z £ 0,2 then Z is T,— equivalent to Z'^ given by ([T|). 
We present an outline of proof of the previous lemma in Section [3l 



At this point it seems natural to establish the following conjecture. 
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Conjecture: For any neighborhood W <Z of Z*™', given by ([I]), and 

for any integer k > there exists Z £ W such that the codimension of Z is k. 

So, based on the conjecture, we have to sharper our analysis. In order to 
get low codimension bifurcation we have to impose some generic assumption. 

Consider Zq = {Xo,Yq) G O'' represented by the model ([T|) with the first 
coordinate of X'^ given by 1 and ki = —1. When r = inv we add the extra 
generic assumption Xq./(p) 7^ on the S— fold point. In [16], Theorem 2, 
we can change the set called Q and to conclude that around the invisible 
S— fold point the vector field Xq is expressed by Xq = (1, —x + oix^) where 
oi / 0. We say that the S— fold point of Xq is contractive (respectively, 
expansive) if ai < (respectively ai > 0). 

According to the previous discussion, we will consider Zq"'',Zq** G 0^ 
written in the following forms: 



Note that X^p^ presents an invisible expansive S— fold point in its phase 
portrait and X^^^ presents a visible one. 

The main question is to exhibit the bifurcation diagram of Zq with either 
T = inv or r = vis. 
We obtain that: 

I- There is a canonical imbedding Fq : M^, — )■ ^0 such that Fq(X, /3) = 
Zl g is expressed by: 



where A G (-1,1), P G (-\/3/2, \/3/2), ai{inv) = -1, ai{vis) = 1, 
a2{inv) = 1 and a2{vis) = 0. Moreover, the two— parameter family given 
by ([4]) is transversal to ^2 and its bifurcation diagram is exhibited (see Fig- 
ures [2T] and [28]) . We observe that there are some typical topological types 




(3) 
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nearby Zq that do not appear in the bifurcation diagram of ZJ^^. For ex- 
ample, when r = inv the configurations in Figures [3] and [Hare excluded and 
when T = vis the configuration in Figure [5] also is excluded. 





Figure 3. 



Figure 4. 




(5) 



Figure 5. 

II- We add an auxiliary parameter // in the following way: 

1 

ai{T){x - X) + a2{T){x - A)2 



Z 




(2/ + /3) 



where A G (-1, 1), (3 
a2{inv) = 1, a2{vis) ■ 



(-V3/2,V3/2), aiiinv) = -1, 
and fi £ {eq, 1) with Eq < 0. 



if y > 0, 



if y < 0, 



ai{vis) = 1, 
By means of 



this latter unfolding its bifurcation diagram cover all topological types near 

^0,0,0- 

The configuration illustrated in Figure [3] plays a very important role in our 
analysis. In this resonant configuration we note, for example, a fold— fold 
singularity or even a loop passing through the saddle point. Only the bi- 
furcation of these two unstable configurations already represents a relevant 
development (the fold— fold singularity was studied recently in [9] and the 
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non— smooth loop bifurcation, as long as we known, was not studied until the 
present work). In fact, this configuration is reached in ([5]), taking fi = fio 
where 

(6) fio = 2- (12/3/(-3 + 6/3 + V9 - 12/32)) 

and A = Ao = (-3 + v^9 - 12/32)/6. 

In this paper we consider just the cases described in Equations ([2]) and 
([3]) , where the first coordinate of X is equal to 1 . When the first coordinate 
of X is equal to —1 a similar approach can be done. 

It is worth mentioning that we detect branches of canard cycles^^ in the 
bifurcation diagram of Z^^^^ p. A canard cycle is a closed path composed by 

pieces of orbits oi X, Y and (see Figures El [8] and [9]). In Section 2 a 
precise definition will be given. 

1.2. Statement of the Main Results. Our results are now stated. Theo- 
rems 1, 2 and 3 are intermediate steps towards Theorem A and Theorems 4, 
5 and 6 are intermediate steps towards Theorem B. Here we follow Definition 
[T] to say when two non— smooth vector fields represent a same topological 
behavior. 

Theorem 1. Take r = inv and fi = fiQ in Equation ([5]), where fiQ is given 
by dl]). Its bifurcation diagram in the {X, /3)— plane contains essentially 19 
distinct topological behaviors (see Fiaure[W\). 

It is easy to see that the cases covered by Theorem 1 do not represent the 
full unfolding of the (Invisible) Resonant Fold— Saddle singularity. Because 
of this, the next two theorems are necessary. Each one of them describes a 
distinct generic codimension two singularity. 

Theorem 2. Take t = inv and /io < /U < 1 in Equation ([5]). Its bifurca- 
tion diagram in the {X, P)— plane contains essentially 21 distinct topological 
behaviors (see Fiaure [2l\} . 

Theorem 3. Take t = inv and Eq < fj, < fio in Equation ([5]). Its bifurca- 
tion diagram in the {X, (3)— plane contains essentially 21 distinct topological 
behaviors (see Figure [2l\) . 

Theorem 4. Take t = vis in Equation (j4]) or equivalently, take r = vis and 
/i = in Equation ([5]). Its bifurcation diagram in the {X, P)— plane contains 
essentially 13 distinct topological behaviors (see Figure \2^) . 

The cases covered by Theorem 4 do not represent the full unfolding of 
the (Visible) Resonant Fold— Saddle singularity. Because of this, the next 
two theorems are necessary. Each one of them describes a distinct generic 
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codimension two singularity. 

Theorem 5. Take r = vis and < fi < 1 in Equation ([5]). Its bifurca- 
tion diagram in the {X, P)— plane contains essentially 13 distinct topological 
behaviors on (see Figure\2B^. 

Theorem 6. Take r = vis and eo < A* < in Equation ([5]). Its bifurca- 
tion diagram in the {\, (3)— plane contains essentially 13 distinct topological 
behaviors (see Figure\2B(j- 

Finally, we are able to state the main results of the paper. 

Theorem A. Equation ([5]) with r = inv generically unfolds the (Invisi- 
ble) Resonant Fold— Saddle singularity. Moreover, its bifurcation diagram 
exhibits 61 distinct cases representing 25 distinct topological behaviors (see 
Figure [2^) . 

Theorem B. Equation ([5]) with t = vis generically unfolds the (Visible) 
Resonant Fold— Saddle singularity. Moreover, its bifurcation diagram ex- 
hibits 39 distinct topological behaviors (see Figure\3^). 

The paper is organized as follows: in Section 2 we give the basic theory 
about Non— Smooth Vector Fields on the Plane, in Section 3 we prove The- 
orem 1, in Section 4 we prove Theorem 2, in Section 5 we prove Theorem 
3, in Section 6 we prove Theorem A and present the Bifurcation Diagram 
of 

Section 7 we prove Theorem 4, in Section 8 we prove Theorem 
5, in Section 9 we prove Theorem 6 and in Section 10 we prove Theorem B 
and present the Bifurcation Diagram of • 

2. Preliminaries 

We distinguish the following regions on the discontinuity set S : 

(i) Si C S is the sewing region if {X.f){Y.f) > on Si . 

(ii) S2 C S is the escaping region if (X.f) > and (Y.f) < on S2. 

(iii) S3 C E is the sliding region if (X.f) < and (Y.f) > on S3. 

Consider Z £ . The sliding vector field associated to Z is the vector 
field Z^ tangent to S3 and defined at g € S3 by Z^{q) = m — q with m being 
the point where the segment joining q + X{q) and q + Y{q) is tangent to S3 
(see Figured]). It is clear that if (7 € S3 then q £ T,2 for —Z and then we 
can define the escaping vector field on S2 associated to Z by Z*^ = — (— Z)*. 
In what follows we use the notation Z^ for both cases. 

We say that (7 G S is a T,— regular point if 
(i) iX.fiq)){Y.f{q)) > or 
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Figure 6. Filippov's convention. 



(ii) {X.f{q)){Y.f{q)) < and Z^{q) 7^ (that is g G E2 U E3 and it is 
not an equilibrium point of Z^). 

The points of S which are not S— regular are called Ti — singular. We dis- 
tinguish two subsets in the set of S— singular points: S* and S^. Any q € Y^p 
is called a pseudo equilibrium of Z and it is characterized by Z^{q) = 0. Any 
g G E* is called a tangential singularity and is characterized by Z^{q) 7^ 
and X.f{q)Y.f{q) = (g' is a contact point of Z^). 

A pseudo equilibrium g' G is a Ti— saddle provided one of the following 
condition is satisfied: (i) g G S2 and q is an attractor for Z^ or (ii) g G S3 
and g is a repeller for Z^. A pseudo equilibrium g G is a Ti—repeller 
(resp. T,— attractor) provided q G S2 (resp. q G S3) and q is a repeller 
(resp. attractor) equilibrium point for Z^. 

Definition 2. Consider Z G fi^ . 

(1) A curve T is a canard cycle ifV is closed and 

• r contains orbit— arcs of at least two of the vector fields X\y,^, 

and Z^ or is composed of a single arc of Z^; 

• the transition between orbit— arcs of X and orbit— arcs ofY hap- 
pens in sewing points; 

• the transition between orbit— arcs of X (or Y ) and orbit— arcs 
of Z^ happens through T,—fold points or regular points in the 
escape or sliding arc, respecting the orientation. Moreover if 
r 7^ S then there exists at least one visible Ti—fold point on 
each connected component ofVOT,. 

(2) Let T be a canard cycle of Z. We say that 

• T is a canard cycle of kind I if T meets E just in sewing 
points; 

• T is a canard cycle of kind II if T = Y^; 

• T is a canard cycle of kind III if T contains at least one 
visible Ti—fold point of Z . 

In Figures and [21 arise canard cycles of kind I, II and III 
respectively. 
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(3) Let T be a canard cycle. We say that T is hyperbolic if 

• T is of kind I and ri'{p) ^ 1, where rj is the first return map 
defined on a segment T with p G T iti 7; 

• T is of kind II; 

• T is of kind III, S^n S^n T = and either T n S C Si U E2 U S* 

orrns c Si ussus*. 







Figure 

7. Canard 
cycle of 
kind I. 



Figure 

8. Canard 
cycle of 
kind II. 



Figure 

9. Canard 
cycle of kind 
III. 



Remark 1. The expression "canard" is used here because these orbits are 
limit periodic sets of singular perturbation problems (see [1] and [6]J. 

Definition 3. Consider Z £ . A closed path A is a Ti— graph if it is 
a union of equilibria, pseudo equilibria, tangential singularities of Z and 
orbit— arcs of Z joining these points in such a way that A n S 7^ 0. As for 
canard cycles, we say that A is a T,— graph of kind 7 i/ A n S C Si, A is 
a T,— graph of kind II if An T, = A and A is a T,— graph of kind III if 

Ans c S2US3. 



In what follows, in order to simplify the calculations, we take ^ = a + 1 
in dni) and obtain the following expression 



(7) 



Xx = 



1 



ai(r)(x - A) + a2(r)(x - A)^ 

2 (y + z^) ^ 



(l+g) ^ ^ i-l+a) , 



if y > 0, 
if y < 0, 



\ 2 ' 2 

where A G (-1,1), /3 G {-V^/2, a G (-1 + eo,0) where Sq < 0, 

T = inv or vis, ai{inv) = —1, ai{vis) = 1, a2{inv) = 1 and a2{vis) = 0. 
When it does not produce confusion, in order to simplify the notation we 
use Z = {X,Y) or Zx,a,f} = {X,Y) instead Z^^^^ = (Xa,^^,/?). 

Since /io is given by ([6]), we obtain that 
(8) ao = 1 - (12/3/(-3 + 6/3 + ^9 - 12/32)). 
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Given Z = {X,Y), we describe some properties of both X = Xx and 

Y = y,,^. 

The real number A measures how the S— fold point d = (^1,^2) = ('^lO) 
of X is translated away from the origin. More specifically, if A < then d 
is translated to the left hand side and if A > then d is translated to the 
right hand side. 

Some calculations show that the curve Y.f = is given hy y = — /3. 

So the points of this curve are equidistant from the separatrices when a = 
— 1. It becomes closer to the stable separatrix of the saddle point S = Sa,j3 = 
(si,S2) when a G (—1,0). It becomes closer to the unstable separatrix 
of S when a £ (—1 + eo,— !)• Moreover, the smooth vector field Y has 
distinct types of contact with S according with the particular deformation 
considered. In this way, we have to consider the following behaviors: 

• Y~ : In this case /3 < 0. So S" is translated to the 7y— direction with 
y > (and S is not visible for Z). It has a visible S— fold point 

e = e„,/3 = (61,62) = ({Sy/3,0) = (61,0) (see Figure [TOD. 

• Y° : In this case /5 = 0. So 5 is not translated (see Figured]). 

• Y"*~ : In this case /3 > 0. So S" is translated to the y— direction 
with y < 0. It has an invisible E— fold point i = ia,i3 = (^1,^2) = 

( (i-a) '^' ^) ■ Moreover, we distinguish two points: h = hp = (hi, /i2) = 
{—(3,0) which is the intersection between the unstable separatrix 
with S and j = jp = (ji,j2) = (/3,0) which is the intersection be- 
tween the stable separatrix with S (see Figure [TT]) . 

In Figure [11] we distinguish the arcs of trajectory ci joining the saddle 
point 5 of y to /i and 02 joining j to the saddle point S of Y . 




3. Proof of Theorem 1 

In [A, B) C S2US3, consider the point C = (Ci, C2), the vectors X{C) = 
{Di,D2) and Y{C) = (^1,^2) (as illustrated in Figure [12]). The straight 
segment passing through C + X{C) and C + Y{C) meets E in a point p{C). 
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We define the C— map 

p: {A,B) S 
z I — > p{z). 

Since S is the x— axis, we have that C = (Ci, 0) and p{C) G M x {0} can be 
identified with points in M. According with this identification, the direction 
function on S is defined by 

H : {A,B) — > R 

z I — > p{z) — z. 



X ^C + Y{C) 




Y 'c + X{C) 



Figure 12. Direction function. 

We obtain that H is a. C^— map and 

• if H[C) < then the orientation of in a small neighborhood of 
C is from B to A; 

• if H{C) = then c G ; 

• if H{C) > then the orientation of in a small neighborhood of 
C is from A to B. 

Simple calculations show that p{Ci) = (-Oi +c^^)^ ^2 +Ci ) conse- 
quently, 

(9) H{c,) = 

Remark 2. If X.f{p) = and Y.f(p) ^ i/ierz, in a neighborhood Vp 
of p in S, the direction function H has the same signal of Di, where 
^ip) = {Di,D2). In fact, X.f{p) = and Y.f{p) ^ are equivalent to 
Di = and £'2 7^ in (l9l) . So, lim Hipi) = Di, where /cn 7^ and 

(D2,E2)^(0,fco) 

P = {PI,V2)- 



Outline of Proof of Lemma 1. Here we construct a S— preserving homeo- 
morphism h that sends orbits of Z = {X, Y) G ^2 to orbits of Z = {X, Y), 
where Z = Z*"'', given by ([1]), the first coordinate of X is equal to 1 and 
ki = —1. The other cases are treated in a similar way. Consider Aq a fixed 
point of the stable separatrix of the saddle point S oi Y (see Figure [TH|) . 
Let Ti be a transversal section of Y passing through Aq. The section Ti 
also is transversal to Y and it crosses the stable separatrix of the saddle 
point 5 of y in the point Bq. Let Ai G Ti be a point on the left of Aq. 
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The trajectory of Y passing through Ai crosses S in a point A2. In the 
same way, the trajectory of Y passing through Ai crosses S in a point B2. 
The trajectory of X that passes through A2 crosses S in a point A^. The 
trajectory of X that passes through B2 crosses S in S3. Consider a fixed 
point of the unstable separatrix of S. Let T2 be a transversal section of Y 
passing through ^44. The section T2 also is transversal to Y and it crosses 
the unstable separatrix of S in the point B^. The trajectory of Y passing 
through crosses T2 in a point A^. In the same way, the trajectory of Y 
passing through i?3 crosses T2 in a point B^. Let A^ S Ti be a point at the 
right of Aq. The trajectory of Y passing through Aq crosses T2 in a point 
Aj. The trajectory of Y passing through Aq crosses T2 in a point Bj. The 
homeomorphism h sends Ti to Ti, the arc of trajectory 71 = Ai Af> to the 
arc of trajectory 71 = Ai B^ and the arc of trajectory 72 = Aq Ay to the 
arc of trajectory 72 = AqBj. Now we can extend continuously h to the 
interior of the region limited by Ti U 71 U T2 U 72. In this way, there exists 
a E— preserving homeomorphism h that sends orbits of Z to orbits of Z. 




□ 

Proof of Theorem 1. In Cases li, 2i and 3i we assume that Y presents the 
behavior Y~ . In Cases 4i, 5i and 61 we assume that Y presents the behavior 
Y^ . In these cases canard cycles do not arise. 

o Case li. di < ei, Case 2i. di = ei and Case 3i. di> ei: The points of 
E outside the interval {di,ei) or {ei,di), according with the case, belong to 
El. The points inside this interval, when it is not degenerate, belong to S3 
in Case li and to E2 in Case 3i. In both cases H{z) > for all 2; S E2 U S3. 
See Figure [HI 

o Case 4i. di < si, Case 5i. di = si and Case 61. di > si: The points of 
S outside the interval (di, si) or (si, di), according with the case, belong to 
El. The points inside this interval, when it is not degenerate, belong to E3 
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A<0 A=0 A>0 




in Case 4i and to S2 in Case 61. In both cases H{z) > for all 2; G S2 U II3. 
See Figure [151 



A<0 A=0 A>0 




Figure 15. Cases 4i, 5i and 61. 



In Cases 7i — 19i we assume that Y presents the behavior Y'^ . Remember 
that a = l-(12/3/(-3+6^+y9 - 12/32)). in what follows we call Lq = 1/12 

(-9-6/3+^9 - 12/32 + ^2^ 15 + ^/9- 12/32 _ 2/3(_3 + 2/3 + ^9 - 12^), 

Ll = -1/2 + v^9- 12/32/6 and L2 = (-9 + 6/3 + v^9 - 12/32 _^ 

+\/2 15 + V9 - 12/32 _^ 2/3(-3 + 2/3 + v^9 - 12/3^)/12. Observe that when 
\ = Lq there exists an orbit— arc of X connecting h and i. When A = Li 
there exists an orbit— arc of X connecting h and j. When A = L2 there 
exists an orbit— arc of X connecting i and j. 

o Case 7i. A < — /3, Case 81. A = — /3, Case 9i. —[3 < A < Lqj Case 
lOi. \ = Lq and Case lli. Lq < A < Li; The points of S outside the 
interval {di,ii) belong to Si. The points inside this interval belong to S3. 
The direction function H assumes positive values in a neighborhood of di, 
negative values in a neighborhood of ii and there exists only one value 
P = Px^a.fi such that H{P) = (in fact, using the software Mathematica, 
we obtain explicitly the value of P, but its expression is too large, so it will 
be omitted). So, by ([9]), the S— attractor P = {P,0), nearby (0,0), is the 
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unique pseudo equilibrium of Z. In these cases canard cycles do not arise. 
See Figure [TBI 

A < -/3 A = -/3 -/? < A < Lo A = Lo Lq < A < Li 




Figure 16. Cases 7i - 111. 



o Case 12i. A = L\: Since \ = L\ there is an orbit— arc 7^ of X 
connecting the points h and j. It generates a S— graph F = 7^ U cr2 U S'U fii 
of kind I. Moreover, since a = ao, where is given by ([8]), there exists a 
non generic tangential singularity at the point d = i. So, the points of 
belong to Si. As the S— fold point of X is expansive, a direct calculus shows 
that the First Return Map r] : {h,d) — )• {h,d) has derivative bigger than 1. 
As consequence, F is a repeller for the trajectories inside it, d = i behavior 
itself like an attractor weak focus and canard cycles do not arise. See Figure 

[m 




Figure 17. Cases 12i, 13i and 14i. 
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o Case 13i. Li < A < L3; The meanning of L3 will be given below in this 
case. The points of S outside the interval (zi, di) belong to Si and the points 
inside this interval belong to S2. The direction function H assumes positive 
values in a neighborhood of di, negative values in a neighborhood of ii and 
there exists a unique value P = Px,a,p such that H{P) =0. So P = (P, 0) 
is a S— repeller. When A is a bit bigger than Li, the First Return Map r] 
has two fixed points, i.e., Z has two canard cycles. One of them, called Fi, 
born from the bifurcation of the E— graph F of the previous case and the 
other one, called F2, born from the bifurcation of the non generic tangential 
singularity presented in the previous case. Both of them are canard cycles 
of kind F Using the software Mathematica we obtain that Fi is a hyperbolic 
repeller canard cycle and F2 is a hyperbolic attractor canard cycle. Note 
that, as A increases, Fi becomes smaller and F2 becomes bigger. When A 
assumes the limit value L3, one of them collides to the other. See Figure [TTl 

o Case 14i. A = L3.' The distribution of the connected components of S 
and the behavior of H are the same of Case 13i. Since A = L3, as described 
in the previous case, there exists a non hyperbolic canard cycle F of kind I 
which is an attractor for the trajectories inside it and is a repeller for the 
trajectories outside it. See Figure [TTl 

o Case 15i. L3 < A < L2, Case I61. A = L2, Case 17i. L2 < A < /?, 
Case I81. A = /3 and Case 19i. A > f3: The points of S outside the 
interval belong to Ei and the points inside this interval belong to 

S2. The direction function H assumes positive values in a neighborhood of 
di, negative values in a neighborhood of ii and there exists a unique value 
P such that H{P) = 0. So, by dH), the S-repeller P = (P,0), nearby (0,0), 
is the unique pseudo equilibrium of Z. In these cases canard cycles do not 
arise. See Figure [HI 




Figure 18. Cases 15i - 19i. 
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Figure 19. Bifurcation Diagram of Theorem 1. 



The bifurcation diagram is illustrated in Figure [T9j □ 

Remark 3. In Cases lli and 15i the ST-bifurcations (as described in [9]j 
arise. In fact, note that the trajectory passing through h can make more 
and more turns around P. This fact characterizes a global bifurcation also 
reached in other cases as shown in this paper. 

4. Proof of Theorem 2 

Proof of Theorem 2. In Cases I2, and 82 we assume that Y presents the 
behavior Y~ . In Cases 42, 62 and 62 we assume that Y presents the behavior 
y". In Cases 72 — 2I2 we assume that Y presents the behavior Y^ . 

o Case I2. di < ei, Case 22- di = ei, Case 82. di > ei, Case A2. 
di < si, Case 52- di = si and Case 62. di > si: The analysis of these cases 
are done in a similar way as the cases li, 2i, 81, 4i, 5i and 61. 

In what follows we call Mq = (-3 - 3a(-2 + a + 2(-l + q)/3) + 
+ v /9(-l + a Y - 12(-1 + a)2/32)/(6(-l + af), Mi = -1/2 + 
+v /9 - 12/376 and M2 = ( -3 + 6/3 - 8a(-2 + a + 2^) + 
+V^9(-1 + aY - 12(-1 + Q)2a2^2)/(g(_i _^ q,)2^^ Observe that when A = 
Mq there exists an orbit— arc of X connecting h and i. When A = Mi there 
exists an orbit— arc of X connecting h and j. When A = M2 there exists an 
orbit— arc of X connecting i and j. 

o Case 72. A < -/3, Case 82. A = -/3, Case 92. -/3 < A < Mq, Case 
IO2. A = Mo and Case II2. Mq < A < Mi; Analogous to Cases 7i - lli 
changing Lq by Mq and Li by Mi. 

o Case 122. A = Mi.' The points of S outside the interval {di,ii) belong to 
Si and the points inside this interval belong to S3. The direction function 
H assumes positive values in a neighborhood of di, negative values in a 
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neighborhood of ii (see Remark[2]) and there exists a unique value P = P\,a,ii 
such that H{P) = 0. So P = (-P, 0) is a S— attractor. Since A = Mi, there 
is an orbit— arc 7^" of X connecting the points h and j. It generates a 
S— graph r = U o"2 U 5 U (Ti of kind I. Since a > ao, where ao is given by 
([8|), it is straight forward to show that the First Return Map defined in the 
interval (/ii, di) C S do not has fixed points. By consequence, F is a repeller 
for the trajectories inside it and canard cycles do not arise. See Figure [20l 



X = Ml Ml < A < ii 




Figure 20. Cases 122, 132 and 142. 



o Case 182 . Ml < A < ii .■ The distribution of the connected components 
of S and the behavior of H are the same of Case 122. Since Mi < \ < ii, 
there is an orbit— arc 75"' of X connecting j to a point k = (A;i,0) G S, 
where ki G {hi,di), for negative time. Also there is an orbit— arc of Y 
connecting A; to a point / = (^i,0) G S, where l\ G (ii, ji), for negative time. 
Repeating this argument, we can find an increasing sequence (/cj)^^^- We 
can prove that there is an interval / C (/c, d) such that r]' = {(pY o ^x)' < 1 
on /. As P is a S— attractor, there is an interval J C (fc, d) such that r]' > 1 
on J. Moreover, using the software Mathematica, we can prove that 77 has a 
unique fixed point Q G {k, d). As consequence, by Q passes a repeller canard 
cycle r of kind I. See Figure [20l This canard cycle born from the bifurcation 
of the S— graph present in Case 122. The expression of rj is too large, so the 
general case will be omitted. For the particular case when a = —1, (5 = 1/2 
and A = —1/2 + ll-v/6/60, the application rj is given by 



rj{x) 



3,3 

4 2 



^1 



+ I + § + 



1 



2 ^ lOv^ 



2x 3 + 2 



1 + 

2 ^ lOv^ 



+ 2x 



A straight forward calculus shows that the unique fixed point of this parti- 
cular rj occurs when X = -^29/2/10. 

o Case 142. A = ii.- Every point of S belongs to Si except the point 
d = i. The canard cycle present in the previous case is persistent for this 
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case (remember that this canard cycle born from the bifurcation of the 
S— graph of Case 122. So, it radius do not tends to zero when A tends to 
ii). So the non generic tangential singularity d = i behavior itself like a 
weak attractor focus. See Figure [20l 

o Case 152. ii < A < M3 and Case I62. A = M^: Analogous to Cases 
13i — 14i changing Li by ii and L3 by M3, where M3 is the limit value for 
with Fi collides to F2. 

o Case 172- M3 < A < A'h, Case I82. A = M2, Case 192- Af2 < A < /3, 
Case 2O2. A = /3 and Case 2I2. A > (3: Analogous to Cases 15i — 19i 
changing L2 by M2 and L3 by M3. 




Figure 21. Bifurcation Diagram of Theorems 2 and 3. 



The bifurcation diagram is illustrated in Figure [2TJ □ 

5. Proof of Theorem 3 

Proof of Theorem 3. In Cases I3, 23 and 83 we assume that Y presents the 
behavior . In Cases 43, 53 and 63 we assume that Y presents the behavior 
Y^. In Cases 73 — 2I3 we assume that Y presents the behavior Y~^. 

o Case I3. di < ei, Case 23. di = ei, Case 83. di > ei, Case 4^. 
di < si, Case 5s . di = si and Case 63. di > si: Analogous to Cases li, 
2i, 3i, 4i, 5i and 61. 

In what follows we consider Mq, Mi, M2 and M3 like in the previous 
theorem. 

o Case 73. A < -/3, Case 83. A = -(3, Case 93. -/3 < A < Mq, Case IO3. 
A = Mq and Case II3. Mq < A < ii; Analogous to Cases 72 — II2 changing 
Ml hyh. 

o Case 123. X = ii: Every point of T,/{d} belongs to Si. In a similar way 
as Case 182, we can construct sequences {ki)i^^ and (Zi)ieN- Since d = i we 



20 



C.A. BUZZI, T. DE CARVALHO AND M.A. TEIXEIRA 



have that ki ^ d and li — )■ d. So d is a non generic tangential singularity 
that behavior itself like an attractor. See Figure [22l 



A = ii n < A < Afi A = Ml 




Figure 22. Cases 123, 183 and 143. 



o Case 183. «i < A < Mi: Analogous to Case 182 except that there 
is a change of stability on P = (-P, 0), which is a S— repeller, and on F, 
which is an attractor canard cycle of kind L This canard cycle born from 
the bifurcation of the attractor non generic tangential singularity present in 
Case 123. See Figure EH 

o Case 143. A = Mi: Analogous to Case 122 except that occurs a change 
of stability on P = {P,0), which is a S— repeller. This fact generates a 
bifurcation like Hopf near P and it appears a hyperbolic attractor canard 
cycle Fi, of kind I, between P and the S— graph F2. See Figure [22l 

o Case 153. Ml < A < M^ and Case I63. A = M3; Analogous to Cases 
152 — I62, changing ii by Mi. 

o Case 173. M3 < A < M2, Case I83. A = M2, Case 19-3 . M2 < A < /?, 
Case 2O3. A = /3 and Case 2I2. A > /?.- Analogous to Cases 172 — 2I2. 

The bifurcation diagram is illustrated in Figure [2TJ □ 

6. Proof of Theorem A 

Proof of Theorem A. Since in Equation ^ we can take a in the interval 
(—00,0), from Theorems 1, 2 and 8 we derive that this equation, with r = 
inv, unfolds generically the (Invisible) Resonant Fold— Saddle singularity. 

Observe that the bifurcation diagram contains all the 61 cases described 
in Theorems 1, 2 and 3. But some of them are S— equivalent and the number 
of distinct topological behaviors is 25. Moreover, each topological behavior 
can be represented respectively by the Cases li, 2i, 81, 4i, 5i, 61, 7i, 81, 9i, 
lOi, 111, 12i, I81, 14i, 15i, I61, 17i, I81, 19i, 122, I82, 142, I23, I83 and 
143. 

The full behavior of the three— parameter family of non— smooth vector 
fields expressed by Equation ([7]), with r = inv, is illustrated in Figure 
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[23] where we consider a sphere around the point (A, /i, (3) = (0, 0, 0) with 
a small radius and so we make a stereographic projection defined on the 
entire sphere, except the south pole. Still in relation with this figure, the 
pictured numbers correspond to the occurrence of the cases described in the 
previous theorems. As expected, the cases 5i and 52 are not represented in 
this figure because they are, respectively, the center and the south pole of 
the sphere. □ 




/3 < < \ >/3 > 

(3 = 

Figure 23. Bifurcation diagram of the (Invisible) Fold— Saddle singularity. 



7. Proof of Theorem 4 

Proof of Theorem 4- Since X has a unique S— fold point which is visible we 
conclude that canard cycles do not arise. 

In Cases I4, 24 and 34 we assume that Y presents the behavior Y~ . In 
Cases 44, 54 and 64 we assume that Y presents the behavior Y^. In these 
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cases, when it is well defined, the direction function H assumes positive 
values. 

o Case I4. di < ei: The points of E inside the interval {di,ei) belong to 
Si. The points on the left of di belong to S3 and the points on the right of 
ei belong to S2. See Figure [241 

o Case 24. di = ei: Here Si = 0. The vector fields X and Y are 
linearly dependent on di = ei which is a tangential singularity. Moreover, 
it is an attractor for the trajectories of Z crossing S3 and a repeller for the 
trajectories of Z crossing S2. See Figure [24l 

o Case 84. di > e\: The points of S inside the interval (ei,di) belong to 
Si. The points on the left of e\ belong to S3 and the points on the right of 
d\ belong to S2. See Figure [JH 




o Case 44. d\ < si: The points of S inside the interval {di, si) belong to 
Si. The points on the left of di belong to S3 and the points on the right of 
si belong to S2. See Figure [25l 

o Case 54. di = si: Here Si = and S is an attractor for the trajectories 
of Z crossing S3 and it is a repeller for the trajectories of Z crossing S2. 
See Figure [25l 

o Case 64. di > si: The points of S inside the interval {di, si) belong to 
Si. The points on the left of si belong to S3 and the points on the right of 
di belong to S2. See Figure [25l 

In Cases 74 — 184 we assume that Y presents the behavior Y'^ . 

o Case 74. di < hi, Case 84. di = hi and Case 9^. hi < di < ii: 
The points of S inside the interval {di,ii) belong to Si. The points on the 
left of di belong to S3 and the points on the right of ii belong to S2. The 
direction function H assumes positive values on S3 and negative values in 
a neighborhood of ii. Moreover, H{/3X/{—l + = and the S— repeller 
P = (/3A/(— 1 + /3),0) is the unique pseudo equilibrium. See Figure [26l 

o Case IO4. di = ii: Here Si = 0. The vector fields X and Y are 
linearly dependent on the tangential singularity di = ii. A straightforward 
calculation shows that H{z) = (1 — /3)/2 7^ for all z G S/jd}. So di = ii 
is an attractor for the trajectories of Z crossing S3 and a repeller for the 
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A < 



A = 



44 

V 




/ 


\ 



54 




A > 



64 




Figure 25. Cases 44, 54 and 64. 

trajectories of Z crossing S2- Moreover, A = {d} U dj U (72 U {S"} U cri U /id 
is a S— graph of kind III in such a way that each Q in its interior belongs 
to another S— graph of kind III passing through d. See Figure [26l 



A < -/3 




104 









Figure 26. Cases 74 - IO4. 



o Case II4. ii < di < ji, Case 12^. di = ji and Case 184. ji < di: 
The points of S inside the interval belong to Si. The points on the 

left of ii belong to S3 and the points on the right of di belong to S2. The 
direction function H assumes positive values on S2 and negative values in 
a neighborhood of ii. Moreover, H{(3\/ [—1 + /3)) = and the S— attractor 
P = (/3A/(— 1 + /3), 0) is the unique pseudo equilibrium. See Figure [271 



0<A</3 A = /3 /3<A 




Figure 27. Cases II4 - 134. 
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The bifurcation diagram is illustrated in Figure [25J 

8. Proof of Theorem 5 

Proof of Theorem 5. The direction function H has a root Q 

(10) ^ 

2(a + 1) 

+ V((-l + a)(l - /?) - A(l + a))2 + 4/3(1 + a)(l + a + A(-l + a))). 

Moreover, H assumes positive values on the right of Q and negative values 
on the left of Q. Note that when a— t-— IsoQ— s-— oo under the line 
{y = 0} and it occurs the configurations showed in Theorem 4. 

In Cases I5, 25 and 85 we assume that Y presents the behavior Y~. In 
Cases 45, 55 and 65 we assume that Y presents the behavior 1"'^. In Cases 
7 5 — 185 we assume that Y presents the behavior Y^. 

o Case I5. di < ei, Case 2^. di = ei, Case 85. di > ei, Case 4^. 
di < si, Case 5^. di = si and Case 65. di > si: Analogous to Cases I4, 
24, 84, 44, 54 and 64 respectively, except that here it appears the S— saddle 
Q on the left of d and e or S. See Figure [29l 

o Case 75. di < hi, Case 85. di = hi, Case 95. hi < di < ii: Analogous 
to Cases 74 — 94, except that here the S— saddle Q appears on the left of di 
and ii. Here P = (p, 0) where 

(11) ^ 

2[a + 1) 

- V^((-l + a)(l - /3) - A(l + q))2 + 4/3(1 + a)(l + a + A(-l + q))). 

o Case IO5. di = ii: Analogous to Case IO4, except that here appear the 
S— saddle Q on the left of di = ii. 



□ 



= {q, 0) where 
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A < (l + a)/3/(l -a) A = (1 + a)/3/(l - a) A > (1 + a)/3/(l - a) 




Figure 29. Cases I5, 2,5 and 85. 



o Case II5. ii < di < ji, Case 12^. di = ji and Case 185. ji < di: 
Analogous to Cases II4 — 184, except that here the S— saddle Q appears on 
the left of di and ii. 

The bifurcation diagram is illustrated in Figure [28j □ 

9. Proof of Theorem 6 

Proof of Theorem 6. The direction function H has a root Q = {q, 0) where 
q is given by (|10p . Moreover, H assumes positive values on the left of Q 
and negative values on the right of Q. Note that when a— )•— IsoQ— )-oo 
under the line {y = 0} and the configurations shown in Theorem 4 occur. 

In Cases Ig, 2q and 3q we assume that Y presents the behavior Y~ . In 
Cases Aq, 5q and Gg we assume that Y presents the behavior Y^. In Cases 
7q — 13q we assume that Y presents the behavior Y^. 

o Case 1q. di < ei, Case 2q. di = ei, Case 3q. di > ei, Case Aq. 
di < si, Case Sg. di = si and Case 6q. di > si, Case 7q. di < hi, Case 8q. 
di = hi, Case 9q. hi < di < ii, Case IOq. di = ii, Case 11q. ii < di < ji, 
Case 12q. di = ji and Case 13q. ji < di: Analogous to Cases I5, 25, 85, 
45, 55, 65, 75, 85, 95, IO5, II5, 125 and 185 respectively, except that here the 
S— saddle Q takes place on the right of di, ei, si and ii when these points 
appear. 

The bifurcation diagram is illustrated in Figure [28l □ 

10. Proof of Theorem B 

Proof of Theorem B. Since in Equation ([7]) we can take a in the interval 
(—00, 0) we conclude that Theorems 4, 5 and 6 prove that this equation, with 
r = vis, unfolds generically the (Visible) Resonant Fold— Saddle singularity. 
Its bifurcation diagram contains all distinct topological behaviors described 
in Theorems 4, 5 and 6. So, the number of distinct topological behaviors is 
89. 

The full behavior of the three— parameter family of non— smooth vector 
fields expressed by Equation ([7]), with r = vis, is illustrated in Figure [30] 
where we consider a sphere around the point (A, n, f3) = (0, 0, 0) with a small 
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/3 < < \ >P > 

/? = 

Figure 30. Bifurcation diagram of the (Visible) Fold— Saddle singularity. 



ray and so we make a stereographic projection defined on the entire sphere, 
except the south pole. Still in relation with this figure, the numbers pictured 
correspond to the occurrence of the cases described in the previous theorems. 
As expected, the cases 5^ and 55 are not represented in this figure because 
they are, respectively, the center and the south pole of the sphere. □ 
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